Rigorous solution of the spin-1 quantum Ising model with single-ion anisotropy 
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We solve the spin-1 quantum Ising model with single-ion anisotropy by mapping it onto a se- 
ries of segmented spin-1/2 transverse Ising chains, separated by the S z = states called holes. A 
recursion formula is derived for the partition function to simplify the summation of hole configura- 
tions. This allows the thermodynamic quantities of this model to be rigorously determined in the 
thermodynamic limit. The low temperature behavior is governed by the interplay between the hole 
excitations and the fermionic excitations within each spin-1/2 Ising segment. The quantum critical 
fluctuations around the Ising critical point of the transverse Ising model are strongly suppressed by 
the hole excitations. 
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The quantum lattice models for which both the ground 
state and the finite temperature thermodynamics can 
be exactly solved are rare. However, these models play 
an important role in the study of quantum criticality. 
A typical example is the uniform S = 1/2 transverse 
Ising model (TIM), which sets a paradigm in elucidat- 
ing the nature of both quantum and thermodynamic 
phase transitions 0, 0, H, 0, HI ■ Experimentally, this kind 
of lattice models can be realized in certain electronic 
materials^, 0, [|| as well as in the optical lattice of cold 
atoms or polar molecules foL ITol] ■ However, in many realis- 
tic situations, the local moments are larger than 1/2 and 
exposed to the interaction of single-ion anisotropy gener- 
ated by the crystal fields[ll[. In particular, in an integer 
spin system, a local spin can be in a neutral polarized 
states with S z = 0. It is unclear how these extra degrees 
of freedom of spins can change the nature of quantum 
criticality revealed by the S=l/2 TIM. To address unam- 
biguously this problem, it is desired to find an extended 
but exactly soluble TIM model with higher spins. 

In this Letter, we study an Ising model of S=l with 
the single-ion anisotropy, defined by the following Hamil- 
tonian 



H = - 



3=1 



[JS*S z 3+l + 2D X (S?) 2 + D Z (S*) 2 ] , (1) 



where L is the lattice length. In the classical limit, i.e 
D x = 0, it reduces to the Blume-Capel model[H, Il3| . 
The ground state of this model has been studied by 
a number of authors [TJ, [l5j]. Oitmaa and Brasch first 
pointed out that in the ground state this model is equiv- 
alent to the S = 1/2 TIM [IH and can therefore be ex- 
actly solved. Here we want to show that the thermody- 
namic quantities of this model can be also rigourously 
calculated. To our knowledge, it is the first quantum 
S = 1 spin model whose thermodynamic properties 
can be rigourously studied without invoking the Bethe 
Ansatz. 

In Eq. JT|), D x and D z are the coupling constants of 



the single-ion anisotropy along the x- and z-axes, respec- 
tively. A single-ion anisotropy along the y-axis can be 
added to this Hamiltonian. However, this term is not in- 
dependent since S 2 = 2. It can be absorbed into the D x 
and D z terms. Eq. ([I]) can be also extended to include the 
spin-1/2 magnetic impurities as well as longitude mag- 
netic field. In the discussion below, free boundary con- 
ditions are assumed. It is straightforward to extend the 
results to the system with periodic boundary conditions. 

At each site of the lattice, Sj can take three values, 
Sj = 0,±1. Effectively, one can regard Sj = ±1 as 
the two polarized spin states of a S=l/2 spin operator 
and Sj = as a hole. A remarkable property of the 
Hamiltonian is that at each site (Sj) 2 commutes with H . 
This means that the hole states (Sj = 0) are decoupled 
from the spin polarized states (Sj — ±1). Thus the total 
number of holes is a good quantum number and can be 
used to classify the eigenstates of H . 

The holes in this system act like non-magnetic im- 
purities. They will separate the system into many in- 
dependent segments of interacting S=l/2 spins. In a 
system of p holes, there are at most p + 1 segments of 
S=l/2 spins. If these holes are located at {xi,---x p } 
with 1 < x\ < ■ ■ ■ < x p < L, it is straightforward to 
show that Eq. ([T]) is exactly equivalent to the following 
Hamiltonian up to a dynamic irrelevant constant (setting 
xq = and x p +\ = L + 1) 



P+i 



H({xi,p}) = MM + Pip, - D a 



(2) 



where 



x n -2 



h(ln) = ~ £ J°3o* +1 - E D *°i> ( 3 ) 

<7p (n = x, z) are the Pauli matrices and l n — x n —x n -\ — 
1 is the segment length. 

The above discussion indicates that in order to di- 
agonalize the model, one needs only to diagonalize the 



2 



Hamiltonian for each individual segment defined by Eq. 
©. The eigenfunction of H is a product of all the 
eigenfunctions of the segment Hamiltonians h(l n ). Cor- 
respondingly, the eigenvalue of H is simply given by the 
sum of the eigenvalues of h(l n ). D z — D x plays the role of 
chemical potential of holes. Thus by adjusting the value 
of D z , one can control the number of holes in the ground 
state. 

h(l) defined by Eq. ([3]) is the Hamiltonian of the S=l/2 
Ising model in a transverse field. It can be reduced to a 
model of non-interacting fermions by the Jordan- Wigner 
transformation. By further diagonalizing this fermionic 
model, the excitation spectrum can be obtained. The 
energy dispersion of the fermionic excitation is given by 



e(l) = ±\D x \y/l + A 2 + 2A cos *, 



(4) 



where A = J/D x and k is determined by the secular 
equation: 



sin(7 + l)k 
sin Ik 



= -A. 



(5) 



This equation is reflection symmetric: if k is a solution, 
then — k is also a solution. For |A| < 1, k has 21 real 
roots within the interval (— n, n]. However, for |A| > 1, 
k has 21 — 2 real roots within (— 7T, 7r] and two opposite 
complex roots. The eigenstates corresponding to these 
two complex solutions of k are two localized states, each 
trapped at one end of the chain. 
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FIG. 1: Minimal energy for exciting a fermion quasiparticle 
Ao or a hole Ah when D z — 0. 

The Hamiltonian (Q]) contains two kinds of excitations. 
One is the fermionic excitation within each segment and 
the other is the hole excitations. The fermion excita- 
tion is gapped except at the critical point |A| = 1. The 
minimal excitation gap is given by 



An = 2\J\ 



The minimal hole excitation gap is determined by the 
minimal energy for creating a hole at one end of the spin 
chain and given by 



E (L-1)-E (L) + D Z -D X 



(7) 



(6) 



where Eq(1) is the ground state energy of h(l). Eq. ([7]) 
holds when Ah > 0. In the case A^ < 0, holes will appear 
in the ground state. These holes will condense and break 
the system into many S=l/2 spin segments. In this case, 
the hole excitation becomes gapless. 

Fig. (JTJ) shows the A dependence of the fermion and 
hole excitation gaps A and A^ for D z — 0. The spec- 
trum of h(L) is unchanged when A changes to —A. Thus 
Ao is symmetric under the refection of A. However, the 
hole excitation gap is non-symmetric when A changes to 
—A. When A > 0, A^ is below Ao except in a narrow 
region around the critical point A = 1. At the critical 
point A = 1, Ah = 0.273 J. In the case A^ < A , the low 
energy excitations are dominated by the hole excitations. 

When D z = 0, there is no hole in the ground state. 
In this case, the ground state can be either in a spin 
ordered phase at which all spins are either ferromagnet- 
ically or antiferromagnetically polarized along the z-axis 
depending on the sign of J, or in a quantum disordered 
phase without any long-range spin order. The transition 
between these two phases at |A| = 1 is of the Ising criti- 
cality. 

Now let us consider the thermodynamic properties of 
the model. At first glance, it seems to be extremely 
difficult to calculate rigourously thermodynamic quanti- 
ties of this system since the holes can take exponentially 
many configurations even though h(l) can be analytically 
solved. However, for the system studied here, we find 
that the partition function can be expressed by the fol- 
lowing recursion formula 

L 

Z = Tr cxp {-/3H) = ^2 ofZ^ (L - p) (8) 

where a — exp [[3{D X — D z )\ and 

ZM(1) = ]T z(h)z(h)-.-z(l p+l ) 

li H — lp+i—1 
I 

= ^Z^Hl -h). (9) 

Zi=0 

Here, we define Z^(l) = z(l) = Tr exp [-j3h(l)}, z(0) = 
1, and z(l) = 2 cosh((3D x ). Thus the partition func- 
tion can be evaluated recursively with Eq. (0), starting 
from a no-hole system. Given z(l), the computing time 
needed just increases quadratically instead of exponen- 
tially with the lattice size. Therefore, this has greatly 
simplified the calculation. It allows us to access readily 
the thermodynamic limit by evaluating exactly the ther- 
modynamic quantities in a sufficiently large lattice, for 
example L = 10 4 , at which the finite size effect can be 
ignored. 



3 



X=1.5 




" '\ *' /' 
V'j'i '•/-1-0 




it I 1 ^ 10 ' ■ 




• A 0.5 / ^^*'~^' ^" -l r; 

^ 



" 6 „ 




(b). 


D 


C(p=0) , 






/ c/^ 






^ dN h /dT 

T^m^m )ILl L t , i i i i i 





2.5 



1.0 
0.5 
0.0 



p. 
H 



0.0 0.4 0.8 



1.2 1.6 2.0 2.4 2.8 
T/J 



FIG. 2: Temperature dependence of (a) the specific heat C 
for several different A and (b) the specific heat C and the 
temperature derivative of the hole excitation number dNh/dT 
for the S — 1 Ising spin model Q with A = 0.5, J = 1 and 
D z = 0. The specific heat for the same model but without 
hole excitations, C(p = 0)), is also shown for comparison. 



T <C min(Ao,Afj), the specific heat drops to zero ex- 
ponentially with decreasing temperature, except at the 
quantum critical points. Above this exponential tem- 
perature dependent regime, two low temperature peaks 
appear in the specific heat curves when < A < 1.5. 
These two peaks result from the hole and fermionic ex- 
citation and appear roughly at T ~ and T ~ Ao, 
respectively. This in fact can be more clearly seen from 
Fig. ((2]-b) where the temperature dependence of the spe- 
cific heat for A = 0.5 is shown and compared with the cor- 
responding curve without any hole excitation C(j> = 0) 
and with the temperature derivative of the thermal av- 
erage number of holes dN h /dT. For A = 0.5, A h < A , 
the low-lying excitations are dominated by the hole exci- 
tations. Thus the low temperature peak arises from the 
hole excitations. 
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From the partition function and its temperature 
derivatives, one can evaluate the free energy and all other 
thermodynamic quantities. However, one can also cal- 
culate directly the internal energy, the correlation func- 
tions, and other measurable variables from the corre- 
sponding segment quantities using the recursion formula 
of the partition function. For example, the internal en- 
ergy is given by 



U = ]TaP(p+l)£ 

p=0 1=0 

+ {D z -D x )N h , 



i(l)z(l)Z( p -V(L-p-l) 



Z{L) 



(10) 



where, Z^Q) = 8 i>Q , u(l) = z' 1 (l)Trh(l) exp[-/3h(l)] is 
the internal energy of a segment, and Nh is the thermal 
average of the hole number defined by 



^ > p=0 



(11) 



The specific heat can then be determined from the first 
order derivative of the internal energy. 

Fig. ||2J| shows the temperature dependence of the 
specific heat for several different A and the tempera- 
ture derivation of the hole excitation number. When 



FIG. 3: The temperature dependence of $(r). The inset 
shows $(T) versus J/T for A = 1 and S = 1 in a wider 
temperature range. 

When D z = 0, there are no holes in the ground state. 
Thus the quantum critical behavior of the model is not af- 
fected by the hole excitations at zero temperature. How- 
ever, at finite temperature, the hole excitations can en- 
hance the thermal fluctuation and suppress strongly the 
onset temperature T* below which the quantum critical- 
ity is observed around the quantum critical point. Re- 
cently, Kopp and Chakravarty calculated the scaled free 
energy coefficient [j| for the S = 1/2 TIM 



$(T) 



2|J|V|A|[F(0)-F(r)] 

rp2 ' 



(12) 



where F(T) is the free energy, and the factor 2\J\ 
is the velocity of the elementary excitations (% = 1) at 
low energy. At the critical point, A = 1, they found 
that <j>(T) already falls into the quantum critical scaling 
regime with $(T) « $(0) = tt/12 at T* < J/2. This sug- 
gests that the quantum criticality in the 5=1/2 TIM 
can persist up to a surprisingly high temperature. How- 
ever, in the presence of holes, we find that the persistence 
of the quantum critical behavior of the system is modi- 
fied. Fig. ([3]) shows <E> as a function of J/T for the S = 1 
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FIG. 4: (a) The phase diagram of the ground state, (b) and 
(c) show the temperature dependence of the magnetization 

m = \J^iji^F^j) I L an d tne corresponding hole excitation 
number Nh for several D z with J — 1, X — 1.5, and Ah(D z = 
0) = 0.447974J. 



TIM at the two critical points (A = ±1) and compared 
with the corresponding results for the 5=1/2 TIM. In 
the case A = 1, $ shows a broad peak around T ~ J/2 
and the onset temperature T* of quantum criticality is 
suppressed by nearly one order of magnitude due to the 
hole excitations. However, at A = — 1, the suppression 



is relatively weak since the hole excitation gap is much 
larger than J/2. 

When the hole excitation has a finite energy gap, the 
ground state is magnetically disordered when A| < 1 but 
ordered (either ferromagnetic or antiferromagnetic, de- 
pending on the sign of J) when |A| > 1. However, for 
sufficiently large but negative D z , denned by Eq. ([7]) 
becomes negative. In this case, the static holes will con- 
dense in low temperatures. Fig. (4-a) shows the ground 
state phase diagram. Across the hole condensed phase 
boundary, both the magnetization m and the hole ex- 
citation number Nh, as shown in Fig. (4-b,c), change 
discontinuously at zero temperature. This is a typical 
characteristic of the first order phase transition. 

In summary, we have developed a recursive method 
to solve the thermodynamics of the 5=1 quantum 
Ising chain with single-ion anisotropy. This allows us 
to evaluate rigorously all thermodynamic quantities in 
the thermodynamic limit. The hole excitations affect 
strongly the low temperature behaviors of the system. 
They enhance the thermodynamic fluctuations and re- 
duce strongly the characteristic temperature of quantum 
criticality. The formula derived in this work hold not 
just for the model studied here. With proper extension, 
they can also be applied to study thermodynamic proper- 
ties of quasi-one dimensional antiferromagnets with non- 
magnetic impurities, such as Sr9(Cui-^.Pd^)Chjl6L fl7l|. 
Cm --r-Zn-^GeOa [18|. as well as other physical systems 
whose Hamiltonian can be written as a sum of indepen- 
dent spin segments, separated by non-magnetic impuri- 
ties. 
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